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Abstract 

The first results from the LHC on jets plus missing energy provide powerful new 
data to test SUSY models. Initial theoretical interpretations of these data have 
concentrated on gravity mediation, usually the CMSSM and its variations. In this 
paper we confront a large class of gauge mediation models with these new data. 
More precisely we consider models of pure general gauge mediation (pure GGM) 
and confront them with the recent experimental results of the ATLAS collaboration. 
We use Herwig++ and RIVET, incorporating the full set of experimental cuts, to 
calculate the signal rates and compare them to the data. Although based on only 
SSpb^ 1 of integrated luminosity, we show that these new data probe and exclude a 
portion of previously allowed parameter space of GGM. 

In addition we investigate the viability of standard SUSY benchmark points, in- 
cluding the Snowmass, CMS and ATLAS sets which encompass other mediation 
scenarios such as gravity, anomaly and gaugino mediation. 



1 Introduction 



Recently the CMS and ATLAS collaborations conducted a first series of searches for 
supersymmetry in 7 TeV proton-proton collisions at the LHC looking for squarks and 
gluinos in final states containing jets and missing energy [lffi]. These analyses were based 



on 35pb~ of data taken in 2010. As no excess above the Standard Model expectations 
was observed in these experiments, their results set limits on the sparticle production and 
consequentially provide stringent new constraints on the allowed regions of parameter 
spaces in SUSY models. So far, most of the theoretical analysis of these initial LHC 
searches has concentrated on constraining the CMSSM and related models [5 10 which 
are expected to arise from gravity mediation scenarios. 

In this paper we will assess the impact of these searches for gauge mediated SUSY 
breaking models. The most stringent constraints on the CMSSMQ currently come from 
the ATLAS zero-lepton searches [3]. Here we will concentrate upon these, and apply them 
to a concise class of General Gauge Mediation (GGM) known as pure GGM 11,12 . The 



results of this analysis are summarized in Fig. [2j We will also comment on constraints for 



ordinary gauge mediation in Section 2.2 



Going beyond gauge mediation, we use the LHC data to test the viability of the 



Snowmass sps 13 , the ATLAS SU 14 and the CMS LM [15 benchmark points, which 



represent a large variety of mediation scenarios (cf. Tab.0. For comparison we combine 
and interpret constraints on the CMSSM, pure GGM and the benchmark points in a 
model independent way in terms of physical squark and gluino masses, in Figs. |5j 



2 Gauge Mediation 



Theories with gauge mediated supersymmetry breaking provide a particularly simple and 
compelling set-up for addressing theory and phenomenology beyond the Standard Model, 
see 16 for a review. On the theoretical side, gauge mediation provides an advantage 
compared to other SUSY-breaking mechanisms (such as gravity mediation) due to its 
automatic avoidance of unsuppressed flavour changing interaction^] Over the last few 
years there has been a surge of interest in gauge mediation which has led to a signifi- 
cant extension and generalisation of its original realisation. The GGM framework, first 



this paper we consider gauge mediation models where the next-to-lightest supersymmetric parti- 
cles (NLSP) are stable on collider timescales. Limits on models where the NLSP decays promptly into 
photons and gravitinos are presented in |4|. 

2 On the other hand, gauge mediation aoes not provide any straightforward WIMP dark matter can- 
didates. Contrary to gravity mediation, the lightest neutralino in gauge mediation will always ultimately 
decay to the gravitino, which is the LSP. This rules out neutralino candidates for dark matter in gauge 
mediation, though not the possibility of gravitino dark matter. 
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introduced in [17) , is suitable for unifying quite general models of gauge mediation in a 
model-independent way. It requires that supersymmetry breaking is communicated to the 
Standard Model (MSSM) sector through gauge interactions at the messenger scale. 



2.1 Pure general gauge mediation 

The GGM formulation (or more precisely pure GGM) is based on the requirement that 
SUSY-breaking effects in the MSSM should disappear in the limit of vanishing Standard 
Model gauge couplings. The resulting description does not require precise knowledge of 
any specific underlying models, which can be weakly or strongly coupled, with explicit 
messenger sectors or direct mediation, or any combination of the above. 

A detailed study of the phenomenology of pure GGM models and their parameter 



spaces was presented recently in 11 12 . As alluded to above, in pure GGM we have 
no direct couplings of the SUSY-breaking sector to the Higgs sector, and therefore the 
soft parameter is approximately zero at the messenger scale. From this starting point 
at the high scale M mess a small but viable value of is generated radiatively at the 



electroweak scale 18 , 19 . Electroweak symmetry breaking then determines the values of 



tan/3 and ji. Since B^ is small, tan/3 is generally large (between 20 and 70). 

The main free parameters of pure GGM models are the gaugino and scalar masses as 



well as the messenger scale 20,21 . In models with messenger fields transforming in com- 
plete and unsplit GUT multiplets, there is a single effective scale Aq for the gaugino masses 
and a single scale A^ for the scalars 21p Generating Aq requires both R-symmetry and 



supersymmetry breaking while A^ is affected only by supersymmetry breaking. For this 
reason and As are two a priori distinct scales in GGM. In the simplest scenario, 
ordinary gauge mediation (on which we will comment in the next subsection), one can 
nevertheless identify these two scales, ~ Ag- 
in GGM the soft supersymmetry breaking gaugino masses at the messenger scale M mess 
are given by 

M~ x .(M mess )= k < M ^ Ag (2.1) 

where ki = (5/3,1,1), kiCti (no sum) are equal at the GUT scale and «j are the gauge 
coupling constants. Similarly, the scalar mass squareds are 

m}(M mcss ) = 2 Ciki A| (2.2) 

where the Ci are the quadratic Casimir operators of the gauge groups. 



3 More generally, if one does not require unification, there are six distinct A scales in GGM: Aa, r , A^, 



with r = 1,2,3 in GGM 17 
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In summary, the value of the high scale M mess , together with and A5 appearing in 



Eqs. (2.1)-(2.2) at M me ss characterise a point in the pure GGM parameter space. In this 
sense pure GGM is the gauge mediation analogue of the CMSSM and mSUGRA models 
with Aq and A5 playing a role similar to the parameters mi and mo in those models. 
However these gravity mediated simple realisations are physically quite distinct from the 
gauge mediated pure GGM framework we analyse here. The main differences in gauge 
mediation include: 



the gravitino is always the LSP and the NLSP can be long lived and is not necessarily 
neutral; 

at the high scale the sfermion masses are not identical, in particular the left- and 



right-handed sfermions have different masses as can be seen from Eq. (2.2); 



M mess is a parameter of the model which is typically much lower than the GUT 
scale. 



Outside the confines of ordinary gauge mediation, where the A^ — As, the pure GGM 
parameter space is populated by many models that predict different values of the ratio of 



gaugino to scalar masses, Ac/ As. This parameter space was investigated in Refs. 11, 12 
from which we adopt Figure [TJ This figure shows the allowed parameter space of pure 
GGM in the (Ac A$) plane (before the new LHC constraints are imposed) for a fixed value 
of M mcss = 10 10 GeV. We used a modified version of SoftSUSY [22], which takes = as 
an input and predicts tan (3 using the electroweak symmetry breaking conditions. Direct 



experimental searches from the Tevatron and LEP (see 12 for more detail) exclude the 
black region in Fig. [T] on the left of the Nessie-shaped pure GGM parameter space, not 
surprisingly it effectively cuts off the lower values of gaugino and scalar masses. Other 
boundaries of the parameter space arise from requiring that there are no tachyons and no 
Landau poles, and that SoftSUSY has not encountered convergence problems during the 
RG evolution between the high and the low scales. 

In gauge mediation models with explicit messengers one expects the ratio Ac/ A$ to be 
close to one (dotted line in Fig. [T]), while for direct mediation models the gaugino masses 
are often suppressed relative to the scalar masses [24}|3~l] (region above the dotted line 
in Fig. [T]). It is also possible to achieve values Aq/As > 1 by increasing the "effective 
number of messengers" 32 37] (region below the dotted line in Fig. [I]). 



The discovery potential for pure GGM models during the early stages of the LHC 
with yfs = 7 TeV was addressed in Ref . [l2] . The left panels in Fig. [2] show three scans 
of the parameter space of pure GGM (taken from [12]), at different messenger scales, 
M mess = 10 8 GeV, 10 10 GeV and 10 14 GeV. The regions expected to be most sensitive 
to pp-scattering at y/s = 7 TeV, correspond to relatively light gluinos and/or relatively 
light squarks. In each figure stop mass contours of 500 GeV and 1 TeV are indicated as 
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Figure 1: Pure GGM parameter space for intermediate messenger scales, M^^s — 
10 10 GeV. The dominant constraints excluding various areas around the allowed (white) 
region are indicated as follows: points in the black region violate the pre-LHC direct 
search limits, while yellow area is excluded by the presence of tachyons in the spectrum. 
In the blue region SoftSUSY has not converged and in the green region a coupling reaches 
a Landau pole during RG evolution. Ordinary gauge mediation lives on the dotted line. 



dotted lines, and the 500 GeV and 1 TeV gluino contours are indicated as solid lines^] 
Furthermore, the diagonal dotted red line corresponds to the boundary between neutralino 
and slepton NLSP. The figures also contain the benchmark points introduced in (12). 



We now implement the new experimental constraints on SUSY searches obtained from 
the ATLAS data on final states with jets, missing energy and no leptons. The details 
of our analysis of this data are explained in section |3j Our results are obtained from 
a Monte Carlo simulation of the signal events using Herwig++ 38,39 and RIVET 40 



implementing all the experimental cuts imposed by ATLAS [3] . Combining the four signal 
regions defined by ATLAS |3| we obtain the constraints shown as the red lines in the right 
panels of Fig. [2} One can clearly see that already the relatively small data sample of 
35 pb _1 provides interesting new bounds on models of pure general gauge mediation. 

As one would expect the excluded regions correspond to relatively low gluino and 
squark masses. This will become clearer in Fig. [6] which shows the GGM exclusion region 



*In the M„ 



10 8 GeV scenario the single clotted contour is for 1 TeV stop masses. 
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Figure 2: The left panels show the pure GGM parameter space in terms of Ac, As defined 
in Eqs. and (Q. From top to bottom we have M mess = 10 8 GeV, 10 10 GeV and 

10 14 GeV. Stop mass contours (500 GeV and 1 TeV ) are indicated as dotted lines, and 
the 500 GeV and 1 TeV gluino lines are solid. The NLSP is neutralino above the diagonal 
red line and stau below. The panels on the right show 95% exclusion contours derived 
from the ATLAS search as red lines, and the black lines indicate uncertainties due to 
scale variations in the NLO cross-section. The colour scale for the right panels shows 
the expected number of signal events normalised to the exclusion limit. The benchmark 



points discussed in 23 are shown as a dot (PGMla middle panel, PGMlb bottom panel), 



triangle for PGM2, a star for PGM3 and finally a square for PGM4. 
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directly in terms of the squark and gluino masses (rather than As and Aq). It is interesting 
to note that there nevertheless exists an allowed narrow wedge shaped region at low values 
of Ag, i.e. for relatively low squarks masses. In this region the NLSP is a stau rather than 
a neutralino. When the NLSP is the lightest stau, the ATLAS jets plus missing energy 
search does not constrain the GGM model. As the NLSP is charged the only missing 
energy in the events comes from the production of neutrinos in the SUSY cascade decays. 
We would expect that the stau would either be reconstructed as a muon, provided that 
its time delay reaching the muon chambers is sufficiently short, or the event is rejected 
due to mismeasurement of the missing transverse energy if its interaction in the muon 
chambers is not recorded. To be conservative we do not consider these events in this 
paper. However we emphasize that this is a very interesting region of parameter space 
as it gives a smoking gun for gauge mediation. This is worthy of further study using a 
similar approach to that in |4lj . 

We can compare the excluded regions with the previously identified best-fit regions 
from a global fit to low energy observables performed in 11 . For both M mess = 10 10 GeV 



and 10 14 GeV the ATLAS search rules out some parameter space which was within the 
95% confidence limits of those fits[^} For M mess = 10 10 GeV a small region of the previous 
68% CL has now been ruled out at low A G and moderate A s . The best-fit points, which 
lie deep in the stau NLSP region, are unaffected by the ATLAS results. This differs from 
recent fits of the CMSSM and related models |6j|, where the best-fit points were more 
strongly affected, although they still remain within the 95% CL of the pre-LHC fits. The 
fact that in PGGM the best-fit points lie in the stau NLSP region further motivates a 
dedicated search for stau NLSPs. 



We find that the lowest viable gluino mass which occurs is 380 GeV, a bound which 
is independent of the messenger scale and the identity of the NLSP. The lowest viable 
squark mass depends on the messenger and on the NLSP identity. For neutralino NLSP, 
we find that the lowest allowed squark mass is 735 GeV which happens when the gluino 
mass is approximately 500 GeV. Since we do not expect the ATLAS search to be sensitive 
to the stau NLSP region, the lowest permissible squark mass there is reduced to 490 GeV 
(where the squark mass is defined as the average mass of the first generation squarks). 
This can be read off from Fig. [6j 



2.2 Ordinary gauge mediation 



Models of ordinary gauge mediation live on the Aq = As slice in GGM parameter space. 
Here, we will also investigate the implications of the ATLAS data on these models. 

In contrast to the discussion of the previous section here we will use a more tradi- 
5 For M mess = 10 s GeV such a fit has not yet been performed. 
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tional approach where tan (3 is treated as a free parameter rather than a prediction. This 
is achieved by deviating from the strict definition of gauge mediation and allowing for 
appropriate non-gauge couplings between messenger fields and the Higgs sectoi|^J This 
can then generate an input value for at the high/messenger scale which is traded for 
tan (3 at the electroweak scale. As a result the essential parameter space is again three- 
dimensional, A, M mess and tan/3. In this scenario the bound from the new LHC data for 
long lived NLSPs is fairly insensitive to the messenger mass with A < 72TeV excluded 
by the ATLAS results (for tan /3 = 10 although the bounds are relatively independent 
of this parameter). This limit entirely arises from the event selection C discussed in the 
next section and corresponds to a squark mass 675 GeV, a gluino mass of 590 GeV and a 
NLSP (lightest neutralino) mass of 94 GeV. 



3 Implementation of the ATLAS event selection and 
analysis of the data 

The ATLAS analysis |3] presents the number of observed events which pass four specific 
event selection criteria, together with the expected number of Standard Model events. 
This can be used, together with an estimate of the number of signal events passing the 
experimental cuts, to determine whether a specific model is ruled out at the 95% con- 
fidence level. Alternatively the limits on the cross sections for non-SM processes given 
in [3] can be used. Here we will do the latter as it already includes most of the statistical 
analysis. 

In order to compare the predictions of a particular BSM model with the ATLAS 
results we therefore need to calculate the expected number of signal events passing the 
cuts in each of the four regions (A,B,C,D) defined in Table 1 of Ref. |3). These regions are 
designed to target light qq, heavy qq, gg and gq production respectively, in the CMSSM. 
This is achieved by imposing different selection criteria on the number of jets (> 2 in A 
and B and > 3 in C and D) as well as on the kinematics (i?™ lss , m e g- and m^), see Ref. [3] 
for details. 

Each SUSY model is a point in the MSSM parameter space which is specified by the 
mass spectrum, SUSY couplings and mixing angles at the electroweak scale. All these 
are contained in SLHA files produced by SoftSUSY [22] starting from the high-scale input 
from GGM or any other model. 

Given the complexity of the signal processes the calculation of the number of signal 



events is best achieved using a Monte Carlo event generator, in our case Herwig++ 38 39 



6 Strictly speaking these extra couplings can also generate new contributions to the Higgs soft masses 
i 2 H and tti 2 h but we will ignore these for simplicity. 
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m [GeV] 

Figure 3: 95% confidence level exclusion limit in the (m , mi) plane for tan (5 = 3, A = 
and /i > in the CMSSM. The solid red line is the result using our signal simulations (the 
solid black lines show the effect of varying the factorization and renormalisation scales 
used to calculate the next-to-leading order SUSY production cross sections by a factor of 
| and 2), whereas the dashed red line is the limit obtained by ATLAS in |3|. The colour 
scale shows the expected number of signal events normalised to the exclusion limit. 



to simulate the signal processes for a given SUSY model. The experimental event selection 
can be implemented using the RIVET |40| analysis framework^] to analyse the hadronic 



final state generated by the Monte Carlo simulation, without the need for a simulation of 
the detector response. 

In principle this is sufficient to calculate the number of expected signal events for any 
new physics model. However, for most new physics signals the matrix elements imple- 
mented in general purpose Monte Carlo event generators are only accurate to leading 
order in perturbative QCD. Herwig+-|- was therefore used to simulate three sets of su- 
persymmetric particle production processes for each point in supersymmetric parameter 
space: a) squark and gluino production, b) the production of an electroweak gaugino in 
association with a squark or gluino and c) the production of slepton and electroweak gaug- 
ino pairs. The fraction of events passing the experimental cuts was then used together 



with the next-to-leading order cross section calculated using Prospino 43-46 to obtain 



the number of signal events passing the cuts for each of the four signal regions. 

In order to check that our simulations and implementation of the experimental cuts 
was reliable we checked that we obtained good agreement with the numbers of events 



7 In addition we used the library based on the results of [42] to calculate the 771^2 variable. 
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passing the cuts for the CMSSM supersymmetric parameter points in (m ,mi) plane for 



tan/3 = 3, A = and \i > supplied as supporting material 47 with |3]. The limit 
we obtain for the CMSSM using our simulations is in good agreement with that obtained 
by ATLASg as can be seen in Fig.[3j These results use the NLO cross sections. As a 
cross check we have also computed the bounds using the leading-order cross sections. The 
resulting limits are somewhat lower but show worse agreement with the ATLAS values. 

The number of events and the 95% confidence level exclusion limit for pure GGM 
with M mess = 10 14 GeV is shown in Fig.[| for each of the four ATLAS-defined regions. 
This exclusion limit is obtained from the maximum non-SM cross sections of 1.3, 0.35, 
1.1 and 0.11 pb, for regions A,B,C and D, respectively given in [3]. As can be seen 
from Fig. [4] the strongest limit is given by the C and D event selections with D giving 
the strongest limit at low and C the strongest limit for high values of Ag. This is 
consistent with expectations from the design purpose of these regions, since at high Ag 
the squark masses are high and the SUSY cross-section is dominated by gg production, 
corresponding to ATLAS' region C At lower values of As we are closer to squark-gluino 
mass degeneracy and so qg production dominates, and region D provides the best search 
limits. The combination of parameter space excluded by the various regions is given in 
the right panels of Fig. [2] for a range of messenger scales. 



Beyond gauge mediation: Analysis of full set of 
benchmark points 



A number of benchmark points have been proposed for the study of supersymmetric 
models at high energy colliders. There is a range of motivations for the selection of these 
points including: satisfying the current experimental constraints, providing the correct 
relic neutralino abundence to satisfy cosmological limits and comparing the potentials of 
different experiments. A number of these points, including the most studied spsla, were 
designed to study the potential of measuring supersymmetric masses at the LHC and 
therefore have relatively light mass spectra. As the LHC data now rule out regions of 
the previously allowed parameter space, nearly half of these points are now excluded. In 
Table [T] we show effects of the ATLAS leptons SUSY search on the Snowmass sps points, 
the ATLAS SU points, the CMS LM points and the PGM benchmarks. The columns of 
the table show the mediation scenario, the calculated cross-sections for each of the regions 
A,B,C and D, and finally whether the point is ruled out by the data and by which regions. 
Six of the ten sps points are ruled out by the ATLAS search, including spsla. Also, four 
of the SU benchmarks are now excluded, and eight of the sixteen LM benchmarks. The 



PGM benchmarks proposed in 12 remain allowed by the recent data. Since the proposed 



8 Except in the high too region where there is a significant scale uncertainty for the signal cross section, 
which we have not included. 
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Figure 4: Event rates of pure GGM (M mess = 10 14 GeV) in the four signal regions cor- 
responding to different sets of experimental cuts |3j. The shading gives the number of 
events predicted in our model, after all cuts have been applied. The red line shows the 
exclusion contour. Below and to the right of the red dotted line the stau is the NLSP 
whereas above and to the left the NLSP is the lightest neutralino. 



benchmarks are mostly at low masses, all the excluded points are ruled out at least by 
region D, if not by more regions. 
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Table 1: Status of SUSY benchmark points. For each point the columns labelled A,B,C 
and D give the cross section for each of the signal regions used in the ATLAS analysis pi. 
The last column shows which of the four regions the point is excluded by using the new 
data. In the GMSB scenerio the NLSP was taken to be stable on collider time scales. The 
starred cross sections are computed at leading-order values whereas all the other values 
are NLO. 
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5 Conclusions 



In this paper we have obtained the first constraints on gauge mediation models of SUSY 
breaking with long-lived NLSP from 35pb _1 of LHC data at 7 TeV involving jets and 
missing energy in the final states. In carrying this out we performed an independent Monte 
Carlo simulation of signal events implementing all experimental constraints imposed by 
ATLAS in (3|. 

Our main results are summarized in Fig. [2] and [6] for gauge mediation and in Tab. [T]and 
Figs. [5] and [7] for a standard set of benchmark points including also non-gauge-mediated 
models and the CMSSM. 

In addition to an interpretation of their results in terms of the CMSSM in [3] ATLAS 
also presented their results as a limit on a squark and gluino masses in a simplifed model 
which only contained squarks, gluinos and a massless neutralino. The latter tends to 
increase the search reach by increasing the energy of the visible decay products. 

In Figs. [5]{7]we interpret and combine the ATLAS constraints on different SUSY models 
by displaying them in terms of physical squark and gluino masses. Fig. [5] shows the 
mapping of the CMSSM mo-mx/2 parameter space into the squark-gluino mass space. We 
note that the region rrig < rriq is not accessible in any of the models we have studied 
due to the effects of the gluino masses on the squark masses during the RG evolution^] 
In general, constraints in the squark and gluino plane are similar in all the models we 
have considered (see Figs.[5]{7]), despite significant differences in the mass spectra of the 
other sparticles. Importantly however this can change if the (N)LSP is not the lightest 
neutralino as can be seen in Fig.[6j Moreover this limit is always weaker than that obtained 
in the simplified model due to the non-vanishing neutralino mass. 



9 While close to one in the CMSSM, the slope of the boundary in (irig,mq) is model dependent. 
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gluino mass [GeV] 



Figure 5: This plot shows constraints on the CMSSM for tan (3 = 3, A = and fi > 
mapped into the plane of the physical squark (average of first generation) and gluino 
masses. The kite-shaped area shows the same region of parameter space as in Fig. [3] The 
grey area is still allowed, whereas the white region inside the kite is now excluded by the 
ATLAS measurements [3]. The region below the diagonal vrtg < rriq is not part of the 
CMSSM parameter space due to the influence of the gluino mass on the squark masses 
during the RG evolution. The dashed green line gives the constraints obtained from a 
simplified model (containing only squarks and gluinos and a massless neutralino) in |3j. 
The reduced sensitivity in the CMSSM is mainly due to the non-negligible neutralino 
mass. The labelled points are the benchmark points of Tab. [T} Red points are now 
excluded whereas green points are still viable. 
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Figure 6: Constraints on pure GGM in the plane of the physical squark and gluino masses 
for three different choices of the messenger scale. For each plot the wedge-shaped region 
shows the previously allowed parameter space, the white part of which is now excluded by 
the ATLAS results. As before, the green dashed line gives the constraints for a simplified 
model. The allowed region near the lower boundary of the model space is due to lack of 
missing energy in models with a stau NLSP. 
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Figure 7: Compilation of constraints in terms of physical gluino and squark masses in 
different models of SUSY breaking. In blue we show constraints in a pure GGM model 
with M mess = 10 10 GeV. The white region enclosed in blue lines is now excluded by the 
LHC data [3], the shaded area is still viable. For comparsion we show in grey the allowed 
and excluded regions for the CMSSM (with tan (3 = 3, A = and fi > 0). We also show 
benchmark points from 12 15] (see Tab. Q. The green points are still allowed, and the 
red ones are now excluded. The dashed green line gives the constraints obtained from a 
simplified model as before. 
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